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Abstract. In the scheme of a quantum nondemolition (QND) measurement, an 
observable is measured without perturbing its evolution. In the context of studies 
of decoherence in quantum computing, we examine the 'open' quantum system of 
a two-level atom, or equivalently, a spin-1/2 system, in interaction with quantum 
reservoirs of either oscillators or spins, under the QND condition of the Hamiltonian 
of the system commuting with the system-reservoir interaction. For completeness, we 
also examine the well-known non-QND spin-Bose problem. For all these many-body 
systems, we use the methods of functional integration to work out the propagators. 
The propagators for the QND Hamiltonians are shown to be analogous to the squeezing 
and rotation operators, respectively, for the two kinds of baths considered. Squeezing 
and rotation being both phase space area-preserving canonical transformations, this 
brings out an interesting connection between the energy-preserving QND Hamiltonians 
and the homogeneous linear canonical transformations. 
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1. Introduction 

In the scheme of a quantum nondemohtion (QND) measurement, an observable is measured 
without perturbing its free motion. Such a scheme was originally introduced in the context 
of the detection of gravitational waves [1]. It was to counter the quantum mechanical 
unpredictability that in general would disturb the system being measured. The dynamical 
evolution of a system immediately following a measurement limits the class of observables that 
may be measured repeatedly with arbitrary precision, with the influence of the measurement 
apparatus on the system being confined strictly to the conjugate observables. Observables 
having this feature are called QND or back-action evasion observables O O ll| • In addition 
to its relevance in ultrasensitive measurements, a QND scheme provides a way to prepare 
quantum mechanical states which may otherwise be difficult to create, such as Fock states 
with a specific number of particles. One of the original proposals for a quantum optical 
QND scheme was that involving the Kerr medium [5], which changes its refractive index as 
a function of the number of photons in the 'signal' pump laser. The advent of experimental 
methods for producing Bose-Einstein condensation (BEG) enables us to make progress in the 
matter- wave analogue of the optical QND experiments. In the context of research into BEG, 
QND schemes with atoms are particularly valuable, for instance, in engineering entangled 
states or Schrodinger's cat states. A state preparation with BEG has recently been performed 
in the form of squeezed state creation in an optical lattice [6]. In a different context, it 
has been shown that the accuracy of atomic interferometry can be improved by using QND 
measurements of the atomic populations at the inputs to the interferometer [7]. 

No system of interest, except the entire universe, can be thought of as an isolated system 
- all subsets of the universe are in fact 'open' systems, each surrounded by a larger system 
constituting its environment. The theory of open quantum systems provides a natural route for 
reconciliation of dissipation and decoherence with the process of quantization. In this picture, 
friction or damping comes about by the transfer of energy from the 'small' system (the system 
of interest) to the 'large' environment. The energy, once transferred, disappears into the 
environment and is not given back within any time of physical relevance. Ford, Kac and Mazur 
[8] suggested the first microscopic model describing dissipative effects in which the system was 
assumed to be coupled to a reservoir of an infinite number of harmonic oscillators. Interest in 
quantum dissipation, using the system-environment approach, was intensified by the works of 
Galdeira and Leggett [^, and Zurek among others. The path-integral approach, developed 
by Feynman and Vernon [11] , was used by Galdeira and Leggett [9] , and the reduced dynamics 
of the system of interest was followed taking into account the influence of its environment, 
quantified by the influence functional. In the model of the fluctuating or "Brownian" motion 
of a quantum particle studied by Galdeira and Leggett [9] , the coordinate of the particle was 
coupled linearly to the harmonic oscillator reservoir, and it was also assumed that the system 
and the environment were initially factorized. The treatment of the quantum Brownian motion 
has since been generalized to the physically reasonable initial condition of a mixed state of the 
system and its environment by Hakim and Ambegaokar [12], Smith and Galdeira [13], Grabert, 
Schramm and Ingold [13], and by us for the case of a system in a Stern-Gerlach potential [15j . 
and also for the quantum Brownian motion with nonlinear system-environment couplings [16] . 

An open system Hamiltonian is of the QND type if the Hamiltonian Hs of the system 
commutes with the Hamiltonian Hsr describing the system-reservoir interaction, i.e., Hsr 
is a constant of motion generated by Hs- Interestingly, such a system may still undergo 
decoherence or dephasing without any dissipation of energy \n\ [T8] . 

In this paper, we study such QND 'open system' Hamiltonians of particular interest in the 
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context of decoherence in quantum computing, and obtain the propagators of the composite 
systems explicitly using path integral methods, for two different models of the environment. 
The aim is to shed some light on the problem of QND measurement schemes. Can one draw 
upon any familiar symmetries to connect with the time-evolution operation of these QND 
systems of immense physical importance? 

We take our system to be a two-level atom, or equivalently, a spin-1/2 system. We 
consider two types of environment, describable as baths of either oscillators or spins. One 
cannot in general map a spin-bath to an oscillator-bath (or vice versa) ; they constitute distinct 
universality classes of quantum environment [121. '^^^ ^^^^ ^^^^ oscillator-bath models 
(originated by Feynman and Vernon [llj ) describes delocalized environmental modes. For the 
spin-bath, on the other hand, the finite Hilbert space of each spin makes it appropriate for 
describing the low-energy dynamics of a set of localized environmental modes. A difficulty 
associated with handling path integrals for spins comes from the discrete matrix nature of the 
spin-Hamiltonians. This difficulty is overcome by bosonizing the Hamiltonian by representing 
the spin angular momentum operators in terms of boson operators following Schwinger's theory 
of angular momentum [20]. 

We then use the Bargmann representation [5T] for all the boson operators. The 
Schrodinger representation of quantum states diagonalizes the position operator, expressing 
pure states as wave functions, whereas the Bargmann representation diagonalizes the creation 
operator 5^, and expresses each state vector in the Hilbert state 7i as an entire analytic 
function f{a) of a complex variable a. The association 1-0) — > fioi) can be written 
conveniently in terms of the normalized coherent states |a) which are the right eigenstates 
of the annihilation operator b: 

b\a) = a\a), 

/ 't \ f ^ \ l\2 -'"112, /* \ 

[a \a) = exp I— — 21*^1 "^j' 

giving 

/(a) = e-HV2 

We obtain the explicit propagators for these many-body systems from those of the expanded 
bosonized forms by appropriate projection. 

The propagators for the QND Hamiltonians with an oscillator bath and a spin bath 
are shown to be analogous to the squeezing and rotation operators, respectively, which are 
both phase space area-preserving canonical transformations. This suggests an interesting 
connection between the energy-preserving QND Hamiltonians and the homogeneous linear 
canonical transformations, which would need further systematic probing. 

The plan of the paper is as follows. In section 2 we take up the case of a QND-type of 
open system Hamiltonian where the bath is a bosonic one of harmonic oscillators. In section 
2.1 we consider a case, which is a variant of the previous one, wherein we include an external 
mode in resonance with the atomic transition and obtain its propagator. In section 2.2 we 
discuss the non-QND variant of the Hamiltonian which usually occurs in the literature in 
discussions of the spin-Bose problem [22t [25] . In section 3 we treat the case of a QND-type 
of open system Hamiltonian where the bath is composed of two-level systems or spins. The 
structure of the propagators in the two cases of the oscillator and spin baths is discussed in 
section 4, and in section 5 we present our conclusions. 
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2. Bath of harmonic oscillators 

We first take the case where the system is a two-level atom interacting with a bosonic bath of 
harmonic oscillators with a QND type of coupling. Such a model has been studied |24[ [25l |26] 
in the context of the influence of decoherence in quantum computation. The total system 
evolves under the Hamiltonian, 



Hi = Hs + Hr + HsR 

^-(Jz + E ^^kblhk + 9k(pk + h\)a,. (1) 



fc=i ^ ^ fc=i 

Here Hs, Hr and Hsr stand for the Hamiltonians of the system, reservoir, and system- 
reservoir interaction, respectively. We have made use of the equivalence of a two-level atom and 
a spin-1/2 system, cTz denote the standard Pauli spin matrices and are related to the spin- 
flipping (or atomic raising and lowering) operators 5+ and S-: Ux = Sj^ + S-., az = — 1. 
In ([T|) bj^,bk denote the Bose creation and annihilation operators for the M oscillators of 
frequency uj^. representing the reservoir, gi^. stands for the coupling constant (assumed real) for 
the interaction of the field with the spin. Since [Hs, Hsr] = 0, the Hamiltonian ([1]) is of QND 
type. 

The explicit propagator exp(— ^) for the Hamiltonian ([1]) is obtained by using functional 
integration and bosonization j221l27j. In order to express the spin angular momentum operators 
in terms of boson operators, we employ Schwinger's theory of angular momentum ^20j by 
which any angular momentum can be represented in terms of a pair of boson operators with 
the usual commutation rules. The spin operators cr^ and ax can be written in terms of the 
boson operators a^j, aj^ and o^, ojj, as 

t t 

t t 
ax = al^ap + a^a^. 

In the Bargmann representation [21] the actions of h and 6^ are 

6t/(a) = a7(a), 

^ (2) 

where \a) is the normalized coherent state. The spin operator becomes 

Here the variable [i* is associated with the spin-down state and the variable 7* with the spin-up 
state. 

The bosonized form of the Hamiltonian (1) is 



fc=i 




huj ^ I ^ d \ f ^ d d 



+ — 1^9k\al + — ][^*—-P*—]. (4) 



Here a^, are the Bargmann representations for 6], and ht, respectively. A particular 

k 

solution of the Schrodinger equation for the bosonized Hamiltonian (4) is 

Ui = Uoop*f3' + C/oi/3*7' + UwYP' + Unl*l', (5) 



Functional integral treatment of some quantum nondemolition systems 



where the amphtude Uij are functions of time as well as the coherent state variables associated 
with the boson oscillators, with the initial condition 



Uij{t = 0) = exp I ala'k | 5ij 



(^J =0,1). 



(6) 



The initial state for the expanded propagator associated with the bosonized Hamiltonian (5) 
is 



U{t = 0) = exp jj^^a^^l exp {/?*/?' + 7*7'} • 



(7) 



If the Hamiltonian is in the normal form given by H{a*, g^;^, the associated propagator is 
given as a path integral over coherent state variables as |28] 

U{a*,t;a',0) = y"D{a}expJ ^ a*(r+)a(r) - ^ ^ dr-ff (a*(r+), a(r), r) I .(8) 



0<T<t 



Here J2 a*(T+)a(r) stands for J2 iii the subdivision of the internal [0,t], 

0<T<t j=0 

i.e., where r stands for Tj, r+ stands for the next point rj-|_i in the subdivision. Also, in the 
subdivision scheme, 

* N-l 

I dTH{a*{T+),a{T),T) = H{a*{Tj+i),a{T^),T^)AT^. 
i=o 



Here the path differential in (8) is 

0<T<t 

where the weighted differential is 



l^(r) = — exp ^— |a(r)p^ (i^a(r). 
Using (8), the propagator for the bosonized Hamiltonian (4) is 

ui(a*,/3*,7*,t;a',/3',y,0) = J {cx}U^ {f3}-D' {-f} 



(9) 



(10) 



X exp < Y 

[ 0<T<t 



M 



k=l 



+ /3*(r+)/3(T)+7*(T+)7(r) 

M * 

Y / dTuJkal{T+)akiT) 

7*(r+)7(T)-/3*(T+)/3(r) 



k=l ] 



i— I dr 



k=i 
r(r+)/?(T 



al{T+) + ak{T) 



7*(r+)7(r) 
(11) 
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In Eq. (11) Q; is a vector with components {a,fc}, and D^jo;} = Hfeii D^{afc}. 

Now we introduce a complex auxiliary field /(t) to decouple the interaction term in (11) 

as 

/ . Af * 

e^P E / drgk [aUT+) + afc(r)] [7*(r+)7(r) - P*{t+)(3{t) 

r _A/ * 

= I D2{/}exp dTr{T)gk{al{T+) + ak{T)) 



X exp 



drfir)- (7*(r+)7(T) - /?* (T+)/3(r)) 



.0 



Here we have used the 5-functional identify, [22j 

t 



J D2{x}P[x*(t)]exp I J dTy{T)x{T) \ = P[y{t)], 



where D^{x} is the functional differential 



D2{x} = exp (- I dT\x{T)A n (— ) d^<r) 



(12) 



(13) 



(14) 



and P[x*(t)] is an explicit functional of x* only. Using (12), the bosonized propagator (11) 
can be written as 

ni(a*,/3*,7*,t;a',/3',7',0) = j D2{/}Gi(a*, t; a', 0; [/*]) 

xiV,(/3*,7*,t;/3',7',0;[/]). (15) 

Here Gi stands for the propagator for 

M 



k=l 

Ni is the propagator for 



d 

^k^l^— + f*{t)gkOil + f*{t)gkak 



da 



H 



7 



(16) 



(17) 



These obey the Schrodinger equations ih-^Gi = HGiGi,ih-^Ni = H^iNi with the initial 
conditions 



Gi{t = 0) = exp l^^afcafcj 



Ni{t = 0) = exp {/?*/?' + 7*7'} • 
The propagator Gi is given by 



(18) 



M 



M 



Gi= exp ^aXe-*"'=*-E 



k=l 



k=l 



ialgk / dTf*{T)e 



-iuJk(t-T) 



+ ia'kgk J dTe~'^>=^r{T) 





t T 



+ 9IJ drj dT'e-'^^^^--'^nT)r{T') 




(19) 
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The propagator A^i is given by 

iVi = exp {QooP*P' + Qoif3*l' + Qm*P' + Qii7*7'} 



1=0 



(r,7*)Q 



7 



Here Q{t) is given by 

Q{t) = exp I Y^zt - ^(Tz j drfir) 



(20) 



(21) 



with Qij{0) = 6ij,Q{0) = I. 

Thus the propagator for the bosonized Hamiltonian (4) as given by (15) becomes 



1=0-^ 



(22) 



The propagator for the Hamiltonian (1) is obtained from (22) by taking the / = 1 term in 
the above equation. By making use of the (^-functional identity (13) the amplitudes of the 
propagator for the Hamiltonian (1) are obtained in matrix form as 



Ul 



Uoo Uoi\ 



( M 

exp < ^ a^a^e' 
U=i 

Afe^ 



where 



A 



B 



k=l 



X e 







\2 M 2 



M 



k=\ ^ 



^9k_ , ^ 
2 Wfe 



(23) 



(24) 



(25) 



(26) 



Here we associate the values a* with time t and a' with time t = as is also evident 
from (8). The simple form of the last term on the right-hand side of (23) reveals the QND 
nature of the system-reservoir coupling. Since we are considering the unitary dynamics of 
the complete Hamiltonian ([T]) there is no decoherence, and the propagator (23) does not have 
any off-diagonal terms. In a treatment of the system alone, i.e., an open system analysis of 
Eq. ([T]) after the tracing over the reservoir degrees of freedom, it has been shown [25] that 
the population, i.e., the diagonal elements of the reduced density matrix of the system remain 
constant in time while the off-diagonal elements that are a signature of the quantum coherences 
decay due to decoherence, as expected. 

Note that though the commonly used coordinate-coupling model describing a free particle 
in a bosonic bath, explicitly solved by Hakim and Ambegaokar [12], with 



M 



H 



(27) 
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is seemingly not of the QND type, it can be shown to be unitarily equivalent to a Hamiltonian 
of the QND type as follows: 

,2 M -, M 



M 



where \J\ and \Ji are the unitary operators 



\J\ = exp 



U2 = exp 



lir I Pj 1 2 



M 



(28) 

(29) 
(30) 



A/ 



The above Hamiltonian ([281) is of the QND type with [Hs,Hsr] = [Py2,P ^ ujjqj] = 0. It 



is commonly known as the velocity-coupling model [29]. 



2.1. An external mode in resonance with the atomic transition 

In this subsection we consider a Hamiltonian which is a variant of the one in (1): 



huj 



Ho = — (Jz + hO.a'^a a 



M 



M 



E hujkb\bk + ^ E 3k{bk + K)(Tz. 



k=l 



k=l 



Here 



n = 2e.d*, 



(31) 



(32) 



where d is the dipole transition matrix element and e comes from the field strength of the 
external driving mode Ei^t) such that 

^^(t) = ee-*"* + e*e^'^*. (33) 

Here we have used the form — ^u^, associated with the external mode, instead of the usual 
form —^Cx and (31) is of a QND type. Proceeding as in section 2 and introducing the symbol 
u* for the Bargmann representation of the external mode a^ we have the amplitudes of the 
propagator for (31) in matrix form as 

f M ) 



U2=[rr TT ) = exp <^ 2^ 

U=i 



X exp 



3B2 







where A is as in Eq. (24), 

M 



B2 = Y<Pk {a*k + a'k) + i 



k=l 



Lo — n 



(34) 



(35) 



and (pk is as in Eq. (26). 
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2.2. Non-QND spin-Bose problem 

In this subsection we consider a Hamihonian that is a variant of the spin-Bose problem 
|22l [27l I23j . This addresses a number of problems of importance such as the interaction 
of the electromagnetic field modes with a two- level atom [30J [31]. Another variant of the 
spin-Bose problem has been used for treating problems of phase transitions [32i 133] and also 
to the tunnelling through a barrier in a potential well [31]. Our Hamiltonian is 



k=l 



M 



(36) 



k=l 



This could describe, for example, the interaction of M modes of the electromagnetic field with 
a two-level atom via a dipole interaction. This has a form similar to Eq. (1) except that here 
the system-environment coupling is via ax rather than cr^. This makes the Hamiltonian (|36p 
a non-QND variant of the Hamiltonian (1). We proceed as in Section II with H^^ (17) now 
given by 

huj 

The propagator for if^r^ (|57|) has the same form as Ni (20) but with Q now satisfying the 
equation 



^ ~ ^ dp* 



1_ — 

d(3* dj* 



(37) 



d_ 

dt 



(38) 



with Qij{0) = 5ij, Q{0) = I. This is solved recursively to yield the series solution 

oo 



Q{t) 



n=0 







X exp 



n 



(39) 



Using Eq. ([39]) and proceeding as before, we obtain the amplitudes of the propagator for the 
Hamiltonian (|36p in matrix form as 

fUoo Uoi 
\Uio Un 

( M 



U3 



exp 



.k=l 



oo 
n=0 



J dTn J dr„_i... J dri exp 



cosh^^X^"^) sinh I^X*-"'*) 

(-l)"sinh(x^"^) (-l)"cosh (x^' 



(40) 
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where 



An) 



(2n + 1) - iwfct + (-l)"+^e 



n+l -iuj^t 



1=1 1=1 

n p—1 
p=2 g=l 



(41) 



and 



X 



in) - _ 



M 

u^^^ gk_ 

2 h ^k 



l+l -iuj^{t-Ti) 



24 E(- 



(42) 



This agrees with the results obtained in |22[l27j. The matrix on the right-hand side of Eq. (40) 
contains diagonal as well as off-diagonal terms in contrast to the matrix on the right-hand side 
of Eq. (23) in which only diagonal elements are present. This is due to the non-QND nature 
of the system-bath interaction of the Hamiltonian described by Eq. (j36p whose propagator is 
given by Eq. (40), whereas Eq. (23) is the propagator of the Hamiltonian given by Eq. (1) 
where the system-bath interaction is of the QND type. The simpler form of the structure of 
the propagator (23) compared to the non-QND propagator (40) reflects on the simplification 
in the dynamics due to the QND nature of the coupling. 



3. Bath of spins 

Now we consider the case where the reservoir is composed of spin-half or two-level systems, as 
has been dealt with by Shao and collaborators in the context of QND systems [18j and also 
quantum computation [35], and for a nanomagnet coupled to nuclear and paramagnetic spins 
|19j . The total Hamiltonian is taken as 

H4 = Hs + Hr + HsR 

= -^5*2 -h ^ /iWfcCJ^fe — ^ CfcCTa-fcS'^. (43) 

k=l k=l 

Here we use 5^ for the system and a^k^crxk for the bath. Since [Hs,Hsr] = 0, we have a 
QND Hamiltonian. In the Bargmann representation, we associate the variable f3* with the 
spin-down state and the variable 7* with the spin-up state for the bath variables, and we have 

* d d 



97* df3* 

Similarly, the bosonization of the system variable gives 

_ d „. d 
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where the variable 6* is associated with the spin-down state and the variable ^* with the 
spin-up state. The bosonized form of the Hamiltonian (43) is given by 

+ ^E«(7j5g. + ffla:^j(r^-«-^). (46) 

A particular solution of the Schrodinger equation for the bosonized Hamiltonian (46) is 
obtained by attaching amplitudes to the polynomial parts in the products 



M 



u, = (r + eW + O n iPk + il) iP'k + I'k) ■ 

k=l 



(47) 



The initial state for the expanded propagator associated with the bosonized Hamiltonian (46) 



IS 



M 



U{t = 0) = exp [6*6' + re'} n Wk + 7^7^ • 



(48) 



fe=i 



Using (8), the propagator for the bosonized Hamiltonian (46) is 



M 



xexp< 

I 0<T<t 



«4(r,r,/3*,7*,*;^',e',/3',7',0)= n / D2{0}D2{e}D2{/3,}D2{7fc} 

e*{T+)e{T)+e{T+)i{T) 

+ /3^(r+)/3fc(r)+7fe(T+)7fe(r) 

C{r+)i{T)-e*{T+)e{T\ 

-i j dTUk 7fe(T+)7fc(T)-/3fe(T+)/3fe(T) 

L 

t 



dT 



e*(r+)e(r)-0*(r+)^(T) 



(49) 



On the left-hand side of Eq. (49), (3*, 7* are vectors with components {(3k} and {7fc}, 
respectively. Now we introduce a complex auxiliary field /(r) to decouple the interaction 
term in (49) as 

exp j drck bUr+)Pk{r) + PUr+hk{r)] [Cir+)ar) - 9*ir+)9{T)] 



j i^2{/}exp -i j dTnr)ckil*k{r+)0k{r)+Ptir+hkir)) 
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X exp 



JdTf{r)'^{C{r+Mr)~0*{r+)9{T)) 



Using (50) in (49) the propagator for the bosonized Hamiltonian (46) becomes 

M 

k=l-' 



xN,^ (/3^, 7^, t; 4, 7^,0; [/*]), 



where Mi is the propagator for 



d \ ihio 



2 V d^* 89* 



Hm, = -TT - + -7rf{t) - 9*^ , 



d 



d 



and Ni is the propagator for 



d 



di* 89* 



a 



d 



8(31 "9jt 



Here the propagator Mi is 



1 



p=i 



pi 



i9*,C)Q{ ^> 



where Q is given by 



Q(t) = exp|y5,t-|5, /dT/(T)|, 







with Qij{0) = Sij,Q{0) =1. 
The propagator Ni^ is 

oo 



(«.7iE)<3M?' 

7jt 



where satisfies the equation 
8 

This equation can be solved recursively to give 

oo 



n=0 



with 



Qfc(O) (Q) ^ Qk{n) ^ 0) = 0, 

t Tn T2 

Q''^^\t) = {iiUka,,rJ dTnJ dTn-l...J dn 


7-/+... + (-ir/) 



X exp -iGx^Ck 



\0 Tl 



Functional integral treatment of some quantum nondemolition systems 



13 



Using Eqs. (54), (56) with p = 1, I = 1, respectively, in Eq. (51) and making use of Eqs. (55), 
(58), (59), (60) along with the (5-functional identity (13), the amplitudes of the propagator for 
the Hamiltonian (43) are obtained in matrix form (in the Hilbert space of Hji) as 



k = l 71=0 





cos(G'^'(")) isin(e'=W) 



^ ' ' ^-i)"isin(e'=W) (-i)"cos(e'^(")) ' ' ^^-^^ 



where 



Qk(n) ^ ^S^CkAn, (62) 

A„ = X](-l)^+'2r, + (-l)"i. (63) 
j=i 

Now if we expand the terms containing Sz, i.e., make an expansion in the system space, 
in Eq. (61) we get terms such as 



e^2-^' cosie^n = cos{-cM r ^ ] . (64) 



Here we have used the fact that 
„5.A _ (e^ o 



Similarly, 



= (o e-^)- 



e^¥^^*.sin(e^(")) = zsin(|c,^„) ("^ ) . (66) 

The above equations have only diagonal elements. We can see from the above equations that 
there are 16 amplitudes of the propagator for each mode k of the reservoir out of which only the 
energy-conserving terms are present due to the QND nature of the system-reservoir coupling. 

4. Discussions 

We look closely at the forms of the propagators (23) and (61) of the QND type Hamiltonians 
([T|) and (43), respectively. In the first case with an oscillator bath, Eq. (23) involves the 
matrix 

\ 

e-^ j ' 

where B is given by Eq. (25). This can be used to generate the following transformation in 
phase space: 

p) " ( e--^) (p. 

It can be easily seen from Eq. (j67p that the Jacobian of the transformation is unity and it is a 
phase space area-preserving transformation. The first matrix on the right-hand side of (j67p has 
the form of a 'squeezing' operation [36j, which is an area-preserving (in phase space) canonical 
transformation coming out as an artifact of homogeneous linear canonical transformations |37| . 
In the second case of a spin bath, Eq. (61) involves the matrix 

/ cose'^^") isine'^^") \ 

\ {-lyn sin e^(") (-1)" cos e'^^") ) ' ^ ' 



(67) 
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where O*'^") is given by Eq. (62). For particular n and k, we write e'^^") as G. For n even, 
the above matrix (1681) becomes 



cos i sin \ _ jeo-. 



■ ■ n n 1-^ (69) 
I sm Id cos fc) ' 

Using the Campbell-Baker-Hausdorff identity [38j this matrix can be shown to transform the 
spin vector a = {(Tx,cry,crz) as 

(70) 



















COS 29 








sin 29 








)l 







a^e'®'^-. (71) 



i.e., the abstract spin vector is 'rotated' about the x-axis by an angle 29. For n odd, (j68p 
becomes (again writing 9'^*-"''' for particular n and k as 9) 

cos 9 i sin 9 \ _ / cos 9 i sin 9 ~ 
-i sin 9 — cos 9 / ^ \ i sin 9 cos 9 

Thus the n-odd matrix is related to the n-even matrix by the spin-flipping energy. The above 
matrix transforms the spin vector a as 

10 

azc'^"- ( ay ] e-'^^-az = e*^ ( cos 29 sin 29 1 I fJy 1 • (72) 










CJ 







,0 sin 29 -cos 29, 

It can be easily seen from the right-hand side of the Eq. (|72p that the determinant of the 
transformation of the spin vectors brought about by the n-odd matrix (1711) has the value 
unity. It is well known that the determinant of a rotation matrix is unity [39]. Thus we see 
that the above transformation has the form of a rotation. Specifically, it can be seen that 

' cos 29 sin 29 \ _ / cos 29 sin 29 \ 

sin 29 - cos 29 / ~ V - sin 29 cos 29 / ' ^ ' 



and 



cos29 sin29\^ _ /cos29 -sin29\ 



sin 29 cos 29; V sin 29 cos 29 )' ^^"^^ 

Here T stands for the transpose operation. From the above it is seen that the matrix 
()68p has the form of the operation of 'rotation', which is also a phase space area- 
preserving canonical transformation [36] and comes out as an artifact of homogeneous linear 
canonical transformations [37j. Any element of the group of homogeneous linear canonical 
transformations can be written as a product of a unitary and a positive transformation |40tl41|. 
which in turn can be shown to have unitary representations (in the Fock space) of rotation 
and squeezing operations, respectively [S^. It is interesting that the propagators for the 
Hamiltonians given by Eqs. ([1]) and (43), one involving a two-level system coupled to a bath 
of harmonic oscillators and the other with a bath of two-level systems, are analogous to the 
squeezing and rotation operations, respectively. 



5. Conclusions 



In this paper we have investigated the forms of the propagators of some QND Hamiltonians 
commonly used in the literature, for example, for the study of decoherence in quantum 
computers. We have evaluated the propagators using the functional integral treatment relying 
on coherent state path integration. We have treated the cases of a two-level system interacting 
with a bosonic bath of harmonic oscillators (section 2), and a spin bath of two-level systems 
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(section 3). In each case the system-bath interaction is taken to be of the QND type, i.e., 
the Hamiltonian of the system commutes with the Hamiltonian describing the system-bath 
interaction. We have shown the commonly occuring free-particle coordinate coupling model 
to be unitarily equivalent to the free-particle velocity coupling model which is of the QND 
type. For the variants of the model in section 2, we have examined (a) the case where the 
two-level system in addition to interacting with the bosonic bath of harmonic oscillators is 
also acted upon by an external mode in resonance with the atomic transition (section 2.1), 
and (b) the non-QND spin-Bose problem (section 2.2), which could be used to describe the 
spin-Bose problem of the interaction of a two-level atom with the electromagnetic field modes 
in a cavity via a dipole interaction. 

The evaluation of the exact propagators of these many body systems could, apart from 
their technical relevance, also shed some light onto the structure of QND systems. We 
have found an interesting analogue of the propagators of these many-body Hamiltonians to 
squeezing and to rotation, for the bosonic and spin baths, respectively. Every homogeneous 
linear canonical transformation can be factored into the rotation and squeezing operations 
and these cannot in general be mapped from one to the other - just as one cannot in 
general map a spin bath to an oscillator bath (or vice versa) - but together they span the 
class of homogeneous linear canonical transformations and are 'universal'. Squeezing and 
rotation, being artifacts of homogeneous linear canonical transformations, are both phase- 
space area- preserving transformations, and thus this implies a curious analogy between the 
energy-preserving QND Hamiltonians and the homogeneous linear canonical transformations. 
This insight into the structure of the QND systems would hopefully lead to future studies into 
this domain. 
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